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If general relativity is the correct theory of physics on large scales, then there is a differential 
equation that relates the Ifubble expansion function, inferred from measurements of angular diam- 
eter distance and luminosity distance, to the growth rate of large scale structure. For a dark energy 
fluid without couplings or an unusual sound speed, deviations from this consistency relationship 
could be the signature of modified gravity on cosmological scales. We propose a procedure based 
on this consistency relation in order to distinguish between some dark energy models and modi- 
fied gravity models. The procedure uses different combinations of cosmological observations and 
is able to find inconsistencies when present. As an example, we apply the procedure to a universe 
described by a recently proposed 5-dimensional modified gravity model. We show that this leads to 
an inconsistency within the dark energy parameter space detectable by future experiments. 

PACS numbers: 98.80.Es,98.62.Sb,98.80.-k 
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I. INTRODUCTION. 

Several cosmological observations, e.g. 0,|3iI1j0;IEIE 

0, , have demonstrated that the expansion of the uni- 
verse has entered a phase of acceleration. The cosmic 
acceleration and the questions associated with it consti- 
tute one of the most important and challenging prob- 
lems for several fields of physics (astrophysics, gravita- 
tion, h igh ener gy p hysics, and fundamental physics), e.g. 
II Elllllll El- A dark energy component that 
represents 2/3 of the entire energy content of the uni- 
verse has been proposed to explain cosmic acceleration 
[isl Ir^ IT1 | . This dark energy component can lead 
to repulsive gravity because of its negative equation of 
state. 

Besides dark energy, several models based on modifi- 
cations to the gravity sector have been proposed recently 
as alternatives to explain the cosmic acceleration, e.g. 

EllSllll. 

These two families of models, dark energy and modi- 
fied gravity, are fundamentally different, and a question 
of major importance is to distinguish between the two 
possibilities using cosmological data. In this paper, we 
propose a procedure that uses different pairs of cosmolog- 
ical observations in order to address this question. The 
procedure is able to find inconsistencies in the dark en- 
ergy parameter space due to an underlying modified grav- 
ity model. 

Furthermore, a burning question that comes to mind 
after an equation of state of dark energy has been de- 
termined from cosmological observations is as follows: Is 
this really a dark energy equation of state or just a results 
obtained because one tried to fit a dark energy model on 
the top of some modified gravity model? The procedure 
proposed also addresses this question. 

The outline of the paper is as follows. In section two, 
we introDuce the basic idea of the consistency test. Next, 



in section three, we recall some of the commonly used 
parameterizations of dark energy and also an example 
of modified gravity models. In section four, we discuss 
constraints on the expansion history while in section five 
we discuss constraints on the growth rate of large scale 
structure. Section six describes our approach. In section 
seven, we describe our implementation of the procedure 
using cosmological observations and present our results. 
We discuss our results and conclude in the last section. 



II. GENERAL RELATIVITY CONSISTENCY 
RELATIONSHIP. 



Einstein's equations of General Relativity relate the 
curvature of spacetime to the matter and energy con- 
tent of the universe, thus describing the cosmological dy- 
namics. Cosmic acceleration affects both the expansion 
history of the universe (given by the Hubble function 
H{z) = i^rl ^ , 1^ 1 , where a is the scale factor and z 

is the redshift) and the rate at which clusters of galaxies 
grow (given by the growth rate of large scale structure, 
D — S{a)/S{ai), where 6{a) = is the overdensity). 

Importantly, the expansion function must be consis- 
tent with the growth rate function via Einstein's equa- 
tions. If the cosmic acceleration is not due to a dark 
energy component in Einstein's equations then the pres- 
ence of significant deviations from the consistency rela- 
tion between the expansion and the growth rate will be 
a symptom of a breakdown of General Relativity at cos- 
mological scales and a hint to possible modified gravity 
models at very large scales. 

In this paper, the consistency test is implemented us- 
ing specific combinations of simulated future cosmologi- 
cal observations. 
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III. DARK ENERGY MODELS VERSUS 
MODIFIED GRAVITY. 

A. Dark Energy parameterization 

We assume that dark energy can be described in terms 
of a cosmic fluid, with an equation of state w = P/ p that 
can vary with redshift. We parameterize the variation of 
w{z) using: 

• w{z) = Wo + wiz if z < 1 and w{z) ~ + wi oth- 
erwise, e.g. and also the Taylor expansion 

• w{a) = Wo + Wa{l - a) + Wb{l - af, e.g. [Mill. 

Given w{z)^ the dark energy density as a function of red- 
shift can be described using the dimensionless function 

^, s . N , . X /"^ l^wiz'))dz' 

Q(z) = p,eW/p<ie(0)=exp3 / , ^ y . (1) 

We assume a spatially flat universe and do not include 
massive neutrinos or isocurvature perturbations. For su- 
pernova (SN la) calculations, we use the parameters 

= {1]a,wo,wi,M}. (2) 

For weak gravitational lensing (WL) calculations, we use: 

^ {il™/i^,fJA,wo,u;i,n5,(Tg",Zp,^5,^r}. (3) 

For CMB, we use: 

= {r!,„/j',f^6/i',flA,u;o,wi,n„(T^",r} (4) 

where Vide = Pde{^) / Pcrit{^) is the dark energy density 
fraction, and Pcrit(O) = SHq/SttG; flmh^ is the physical 
matter density; flbh^ is the physical baryon density; t is 
the optical depth to the surface of last scattering; dg " 
is the amplitude of linear fluctuations; ris is the spec- 
tral index of the primordial power spectrum; Zp is the 
characteristic redshift of source galaxies for weak lensing 
[33. Isgf ; and are the lensing absolute and relative 
calibration parameters as defined in . 

B. Modified gravity models 

In this paper, we are interested in using a good example 
of modified gravity model but only in order to illustrate 
the procedure proposed. We are not interested here in 
evaluating any particular model of modified gravity. 

Although several models of modified gravity have been 
proposed recently as alternatives to dark energy, we 
chose in this analysis to use the Dvali-Gabadadze-Porrati 
(DGP) model. The model is inspired by higher dimen- 
sional physics and is consistent with current observations, 
e.g. H^IH^- As mentioned, our interest is to use the 
DGP model as an example of deviation from General 
Relativity in order to test the proposed procedure. We 



refere the reader to studied dedicated to the DGP model 
phenomenology HI HI El El 

We provide in this sub-section a brief introduction to 
the model, but we refer the reader to 0, for a full 
descri ptio n. The action for the five-dimensional theory 
is [H I23 

+ J '^^^V~5(4)-R(4) + Smatter, (5) 

where the subscripts 4 and 5 denote quantities on the 
brane and in the bulk, respectively; Af(5) is the five di- 
mensional reduced Planck mass; Mf^^^ — 2.4 x lO^^GeV 
is the four dimensional effective reduced Planck mass; R 
and g are the Ricci scalar and the determinant of the met- 
ric, respectively. The first and second terms on the right 
hand side describe the bulk and the brane, respectively, 
while Smatter IS the action for matter confined to the 
brane. The two different prefactors and 
in front of the bulk and brane actions give rise to a char- 
acteristic length, scale |2^ 

r, = Mf4)/2Mf5). (6) 

If A/(5) is much less than M(4) , then the brane terms 
in the action above will dominate over the bulk terms 
on scales much smaller than Tc, and gravity will appear 
four dimensional. For example, nonrelativistic small- 
scale gravity will obey the Newtonian inverse-square force 
law. On scales larger than Tc, the full five dimensional 
physics will be recovered, and the gravitational force law 
will revert to its five dimensional 1 /r^ form. This is usu- 
ally discussed in terms of gravity leakage into an extra 
dimension. Ref. [23| shows that tuning M(5) to about 
10 — lOOMeV, implying Vc ~ H^^, is consistent with cos- 
mological data. They have also been discussed in [3ll |. 
Low redshift cosmology in DGP brane worlds was stud- 
ied in |[20L |2^ . Defining the effective energy density 

^'■^ = (32^Gr2)' 
Friedmann's first equation becomes 

^Lp + ^ = — (Vp + P^+Vp^) ■ (8) 

For each perfect fluid i on the brane, we may define 
ili = pifl/ Pcritfi- Following |23|, we also define Q.r^ = 
jr^^Hf^^. We focus here on a flat universe contain- 
ing only nonrelativistic matter, such as baryons and cold 

dark matter, in which flj.^ — {^—j^) ■ In this model, 
the gravitational "leakage" into the fifth dimension, on 
large length scales, becomes a substitute for dark energy. 
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FIG. 1: Supernova Hubble diagrams for several dark en- 
ergy and DGP models. Note that the ACDM model (red 
solid line) and the flm ~ 0.20 DGP model (blue dotted) have 
nearly identical Hubble diagrams, but different growth factors 
as shown in FigH The same is true of the SUGRA (green 
dashed) and Qm = 0.27 DGP (black double dotted) models. 



FIG. 2: Growth factors of linear density perturbations for 
several dark energy and DGP models. Comparisons among 
several dark energy and DGP models growth factors of linear 
density perturbations. Note that the growth factor in the 
f^m ~ 0.27 DGP model is suppressed with respect to that in 
the ACDM model, which has the same Urn- 



IV. CONSTRAINTS FROM THE EXPANSION 
HISTORY USING SUPERNOVAE OF TYPE lA. 

For a spatially flat universe with dark energy, the Hub- 
ble parameter H{z) can be expressed in terms of its cur- 
rent value Hq and the function Q{z) as 



(z) = HaVn„,il + zf + il-n,n)Qiz). (9) 

Since probes such as the SN la surveys give information 
about the redshift variation of H(z) but not its current 
value, we define the dimensionless Hubble parameter for 
dark energy models as 



Hvsiz) = Hiz)/Ho = ^nrn{l + z)3 + {l-ara)Q{z). 

(10) 

For spatially flat DGP model, the expansion is given by 



nvgviz) = i(i-a„) + y^(i-a«)2 + r!™(i + z)3. 

(11) 

To probe the expansion history, we consider calibrated 
type la supernovae as standard candles to measure the 
luminosity distance as a function of the redshift. Recall 
that the SN la apparent magnitude as a function of red- 
shift is given by 



m(z) = Slogio {DLiz))+M. 



(12) 



Here, Ai depends on Hq as well as the absolute magni- 
tude of type la supernovae. We treat as a nuisance 
parameter. Meanwhile, the dimensionless luminosity dis- 
tance is 

Dl{z) = {1 + z) f dz'/n{z'). (13) 
Jo 



As well known, m(z) depends on wq and wi only through 
a complicated integral relation, which "smears out" dark 
energy information . The result is a degeneracy which 
allows very different models to have nearly identical 
m(z). Fig^shows m{z) for a few dark energy and DGP 
cosmologies. The figure shows that some DGP models 
are very degenerate with certain dark energy models, 
such as ACDM models or SUGRA-inspired models. Ev- 
ident in the expression for m(z) is another degeneracy, 
between and the equation of state parameters. We 
break this degeneracy by combining SN la data with the 
CMB power spectrum. We assume the availability of 
a sample of 2000 supernovae, evenly distributed in red- 
shift between ZMin = and ZMax = 1-7, with a mag- 
nitude uncertainty per supernova of cr„i = 0.2. As in 
23, 24], we also include a peculiar velocity uncertainty 
ay = 500km/s, and a systematic uncertainty Sm — 0.02 
in bins of Az = 0.1. 



V. CONSTRAINTS FROM THE GROWTH OF 

LARGE SCALE STRUCTURE AND 
GRAVITATIONAL LENSING TOMOGRAPHY. 



For dark energy models, the suppression of the growth 
of large scale structure is due to an increase in the ratio of 
dark energy density to matter density. We calculate the 
growth factor for these models by numerical integration 
of the differential equation from |4Qi. .41.J , 
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3 u'(a) 



2 1 + X{a) 



^.^1X^^ = 0, (14) 
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FIG. 3: CMB power spectra for several dark energy and DGP 
models: ACDM model is in red solid line; SUGRA model is in 
green dashed line; flm ~ 0.27 DGP model is in black double 
dotted line; flm = 0.20 DGP model is in blue dotted line; 
Following l3i|, a modified version of CMBFAST ^ was used 
to generate the DGP plots. 
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FIG. 4: Lensing convergence power spectra for several dark 
energy and DGP models: ACDM model is in red solid line; 
SUGRA model is in green dashed line; flm ~ 0.27 DGP model 
is in black double dotted line; Qm ~ 0.20 DGP model is in 
blue dotted line; 



where ' = d/da, G — D/a is the normaUzed growth fac- 
tor, 



Xia) 



(15) 



and Q{a) is as given earlier. 

For DGP models, the suppression of the growth is due 
to weakened gravity resulting from "leakage" into the 
extra dimension. Following |3lLl3^ . we use 

6 + 2H^^J-A7TGp{l + j^)S^0 (16) 



where 



(17) 



FigIS shows that, compared to a ACDM model with the 
same matter density, a DGP model has a distinct sup- 
pression of the growth factor. Also, Fig[21 displays how 
the degenerate models of Fig^show distinct growth fac- 
tor functions. In the next two sections, we will take 
advantage of this difference in order to detect possible 
signatures of DGP models. 

We consider weak gravitational lensing (cosmic shear) 
tomography as a probe of the growth factor of large scale 
structure. Indeed, cosmic shear observations capture the 
effect of the cosmic acceleration on both the expansion 
and the growth factor. We follow the formalism and con- 
ventions as used in p^. Is^ l and we use the convergence 
power spectrum as our statistic Q| i 
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Pi — ^^O^m 



XH 



Hx) 

a2(x) 



3D 



I 



sinK(x)' 



X]dx, (18) 



where P^o is the 3D nonlinear power spectrum of the 
matter density fluctuation, i5; a{x) is the scale factor; 
and sin/f(x) = K^^^^ sm{K^^'^x) is the comoving angu- 
lar diameter distance to x (for the spatially flat universe 
used in this analysis, this reduces to x). The weighting 
function ^(x) is the source-averaged distance ratio given 
by 



six) 



X sinK(x) 



X) 



dx' 



(19) 



where n{x{z)) is the source redshift distribution (we use 
the distribution n{z) = ^e"^/^", which peaks at 

Zp = 2zo). We integrate numerically the growth factor 
function given by Eqns. 1141 or [T?)l above, and we use the 
linear-to-nonlinear mapping procedure halofit . We 
used the halofit for the DGP models as well because we 
are only interested into using them as an example of de- 
viation from general relativity in order to test our proce- 
dure. Dedicated studies to DGP phenomenology should 
consider using a better approximation for these models. 
As mentioned above, the expansion history is contained 
in the window function, while the growth factor is con- 
tained in the 31? nonlinear matter power spectrum. Fig^ 
shows convergence power spectra corresponding to dark 
energy models and modified gravity (DGP) models. 

The separation of source galaxies into tomographic 
bins improves significantly the constraints on cosmolog- 
ical parameters, and particularly those of dark energy 
od\ . The constraints obtained from different bins are 
complementary, and add up to reduce the final uncer- 
tainties on the parameters. Based on recently proposed 
surveys (e.g. [43), we use a reference survey with a sky 
coverage of 10%, a source density n = 100 gal/arcmin , 
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1 /2 

an rms intrinsic ellipticity {"ffnt) ~ 0.25, and a me- 
dian redshift Zmed = 1-5. We assume a reasonable pho- 
tometric redshift uncertainty of a{z) = 0.05, and we 
split the sources into 10 tomographic bins over the range 
0.0 < z < 3.0 with intervals of Az — 0.3. We have used 
^max = 3000, since the assumption of a Gaussian shear 
field, underlying the Fisher formalism and the halofit 
approximation, may not be valid for larger £s. For the 
minimum I, we take the fundamental mode approxima- 
tion, i.e., we consider only lensing modes for which at 
least one wavelength can fit inside the survey area. 



VI. THE FUNDAMENTALS OF THE 
CONSISTENCY TEST: CONSTRAINTS FROM 
THE EXPANSION VERSUS CONSTRAINTS 
FROM THE GROWTH OF LSS 

The cosmic acceleration affects cosmology in two ways: 
i) It affects the expansion history of the universe and ii) 
It supresses the growth rate of large scale structure in 
the universe. The idea that we explore is that, for dark 
energy models, these two effects must be consistent one 
with another because they are related by General Rela- 
tivity. The presence of significant inconsistencies between 
the expansion history and the growth rate would be the 
signature of some modified gravity models at cosmolog- 
ical scales. We propose a procedure that detects such 
inconsistencies when they are present. 

The approach we explore is as follows. We assume 
that the true cosmology is described by a DGP model, 
and then ask what contradictions arise when the data are 
instead analyzed based on the assumption of a dark en- 
ergy model (as mentioned earlier, we are not particularly 
interested in the DGP cosmology here but we want to 
use it as an example to illustrate the procedure.) 

Because we will generate the data using the DGP 
model, the consistency relation from General Relativity 
between the expansion history and the growth factor of 
large scale structure will be broken. The dark energy 
equation of state Wexpiz) which best fits measurements 
of the expansion will not be consistent with the equa- 
tion of state Wgrowth{z) which best fits measurements the 
growth. 

The methods and steps we use are as follows. 

i) We simulate the data for the expansion and the growth 
using a fiducial DGP cosmology (for example, supernova 
and weak lensing data); 

ii) We use a minimization method in order to find 
the best fit dark energy model to measurements of the 
expansion (the minimization method was discussed in 
detail in [23 and was shown to give similar results to 
those from Monte Carlo Markov Chain method although 
the x^ minimization can have a better handle on degen- 
eracies psjl: 

Hi) we use the x^ minimization in order to find the best 
fit dark energy model to measurements of the growth of 
structure. 



iv) Next, we use a standard Fisher matrix approach to 
calculate the confidence regions (or x^ contours) around 
the two best fit dark energy models (this standard pro- 
cedure is discussed in detail in references [sijl: 

v) We look for significant differences between the two ef- 
fective dark energy parameter spaces. 

These differences will signal an inconsistency between the 
expansion and the growth factor. The source of the in- 
consistency is from point (i) where the data was gener- 
ated using a DGP model, i.e. from our hypothesis that 
the true cosmology is that of a modified gravity DGP 
model. 

The basic idea explored here was also discussed in ref- 
erences 1^ but here we propose to implement the 
idea using pairs of cosmological data sets. In the next 
section, we will show how a procedure using observa- 
tions from supernovae, weak lensing tomography, and the 
CMB allows one to achieve an observational implemen- 
tation of the consistency test. 

Note that our approach is totally different from that 
of jsJIS^I. Their work defined 16 new parameters (ad- 
ditional to the cosmological parameters) describing the 
distance and growth factor as functions of redshift. A 
weak lensing survey and CMB power spectra were then 
used to constrain all 16 parameters simultaneously in or- 
der to search for inconsistencies. 

Our work requires no new parameters and we explore in- 
consistencies between constraints obtained from different 
pairs of cosmological probes as we illustrate in the next 
section. 



VII. A PROCEDURE TO DETECT SIGNATURE 
OF MODIFIED GRAVITY USING [SN lA+CMB] 
VERSUS [WEAK LENSING+CMB] 
OBSERVATIONS. 

We provide here an implementation of the consistency 
test using simulated cosmological observations. This 
observational test allows one to distinguish between dark 
energy models and modified gravity DGP models. We 
assume the true cosmology to be a DGP model and we 
analyze two combinations of simulated data sets using 
dark energy models. We use the methods and steps 
indicated in the previous section. The procedure we use 
is as follows. 

i) We use a fiducial DGP model (with VLm = 0.27) and 
generate supernova magnitudes, weak lensing conver- 
gence power spectrum, and CMB power spectrum. 

ii) Then, we determine the dark energy model that is 
the best fit to the supernova magnitudes and to the 
CMB spectrum generated using the fiducial DGP model. 
Hi) Next, we determine the dark energy model that is 
the best fit to the fiducial weak lensing spectrum and to 
the CMB spectrum. 

iv) Then, we compare the allowed regions in the 
{Q,de,WQ,wi} parameter space between the two data 
combinations in order to look for inconsistencies. 



6 




FIG. 5: Equations of state found using two different com- 
binations of data sets. Solid contours are for fits to SN la 
and CMB data, while dashed contours are for fits to weak 
lensing and CMB data. The significant difference (inconsis- 
tency) between the equations of state found using these two 
combinations is a signature of the DGP model and should 
be detectable by future experiments described in sections IV 
and V. The inconsistency is an observational detection of the 
underlying modified gravity DGP model (assumed here to 
generate the data). 



FIG. 6: Dark Energy parameter spaces {fide, wq, wi) found 
using two different combinations of data sets. The ellipsoid 
to the left is for fit to SN la and CMB data, while ellipsoid 
to the right is for fits to weak lensing and CMB data. The 
significant difference (inconsistency) between the parameter 
spaces found using these two combinations is a signature of 
the DGP model and should be detectable by future experi- 
ments described in sections IV and V. The inconsistency is 
an observational detection of the underlying modified gravity 
DGP model (assumed here to generate the data). 



The supernova combination contains information on 
the effect of the acceleration on the expansion history, 
while the weak lensing combination contains information 
on the growth factor in addition to the expansion. This 
translates into a difference between the two allowed re- 
gions in the dark energy parameter space. 

Our results in Figs. El El and [3 show that the two 
allowed regions in the dark energy parameter space are 
significantly different. This signals the expected incon- 
sistency between the expansion history and the growth 
rate of large scale structure. 

VIII. DISCUSSION 

The source of the inconsistency in Figs. 00 and His 
that we generated the data using a fiducial DGP model 
but then we used dark energy models to fit the data. The 
inconsistency is a consequence of the fact that the fidu- 
cial cosmological model is DGP, and therefore it consti- 
tutes an observational signature of the underlying mod- 
ified gravity model. Fig[7| shows that the inconsistency 
between the two equations of state persists even when 
we consider a third term in the Taylor expansion of the 
equation of state (section IIIA) . This indicates that the 
test is robust to the functional form used for the equation 
of state. 

However, we stress that we did not consider in this 
study dark energy models with couplings or with an un- 



usual sound speed. The study of the effect of the sound 
speed of dark energy on the procedure and also the inclu- 
sion of other models of dark energy and modified gravity 
will be considered in subsequent work. 

Also, another point worth exploring in future work is 
the effect of systematic uncertainties. If such a proce- 
dure is performed using real data and an inconsistency is 
found, then one has to develop methods to test whether 
the inconsistency is due to physics or to systematics in 
either of the data sets. Further work will be necessary in 
order to make this type of test robust and generic. 

In summary, the comparison of measurements of ex- 
pansion history and measurements of the growth rate of 
structure tests the behavior of gravity on large scales. We 
proposed a procedure that uses different pairs of cosmo- 
logical data sets in order to explore this comparison and 
to distinguish between some models of dark energy and 
modified gravity as the cause of the cosmic acceleration. 
Being able to distinguish between the two possibilities 
is an important step in the quest to understand cosmic 
acceleration. 
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